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We study the response of the QCD vacuum to an external abelian chromomagnetic field in the 
framework of a non local Nambu-Jona Lasinio model with the Polyakov loop. We use the Lattice 
results on the deconfinement temperature of the pure gauge theory to compute the same quantity 
in the presence of dynamical quarks. We find a linear relationship between the deconfinement 
temperature with quarks and the squared root of the applied field strength, gH , in qualitative (and 
to some extent also quantitative) agreement with existing Lattice calculations. On the other hand, 
we find a discrepancy on the approximate chiral symmetry restoration: while Lattice results suggest 
the deconfinement and the chiral restoration remain linked even at non-zero value of gH, our results 
are consistent with a scenario in which the two transitions are separated as gH is increased. 



I. INTRODUCTION 

Quantum Chromodynamics (QCD) is nowadays re- 
garded as the theory of strong interactions. The study of 
the QCD vacuum and of the QCD phase diagram is one 
of the most intriguing research topic in modern physics. 
A deep understanding of the vacuum would allow us to 
encompass the mechanisms on which color confinement 
and spontaneous chiral symmetry breaking lie. 

The major knowledge on the QCD phase diagram 
arises from Lattice QCD at finite temperature and zero 
(or small) baryon chemical potential, see for exam- 
ple [U, [3, 3 and references therein for recent results. 
At finite chemical potential Montecarlo simulations with 
three colors suffer the well known sign problem Q . In or- 
der to circumvent this problem several approaches to Lat- 
tice calculations have been suggested 0,11, 0, H, S O] ■ 
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Beside them, the Nambu-Jona Lasinio model 
a very popular tool for making predictions on the QCD 
phase diagram in regions unaccessible by Lattice QCD. 

Spontaneous chiral symmetry breaking is characterized 
in terms of the non vanishing and non perturbative chi- 
ral condensate (ipip) , which is an order parameter in the 
massless quark case. On the other hand, the color con- 
finement is well defined in terms of an order parameter 
only in the pure gauge theory. In this case it is well 
known that the Polyakov loop L can be regarded as the 
order parameter for the confinement-deconfinement tran- 
sition As a matter of fact, L is a gauge singlet field 
that transforms non trivially under a transformation. 
Thus in the phase with (L) =/= the symmetry is spon- 
taneously broken, while it is not broken if (L) = 0. Since 
(L) = in the confined phase and (L) 7^ in the decon- 
fined system, one argues that deconfinement-confinement 
transition can be described as a restoration of the spon- 
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tancously broken Z3 symmetry. Dynamical quarks break 
the Z3 symmetry, thus the Polyakov loop in presence of 
dynamical quarks is no longer a good order parameter for 
deconfinement. Nevertheless lattice measurements show 
the existence of a Z3 crossover in the range of tempera- 
ture in which both chiral condensate and particle suscep- 
tibility have a crossover too, denoting that the Polyakov 
loop is still a good indicator of the deconfinement tran- 
sition. 

In this paper, we use an external abelian chromomag- 
netic field to probe the vacuum structure of QCD. This 
is a trick that allows us to learn something more about 
the QCD vacuum, by studying its response to external 
perturbations. This approach has revealed to be success- 
ful in Lattice QCD calculations [3, [H, ■ One of the 
main results obtained is that the deconfinement temper- 
ature at zero baryon density decreases as the strength of 
the applied field gH is increased. In more detail, in the 
case of the pure glue theory, Lattice data on the decon- 
finement temperature Tc as a function of y/gH are well 
interpolated by a linear fit for gH lower than a critical 
value gHc, and consistent with zero for gH > gHc- In 
presence of dynamical quarks. Lattice data are consistent 
with a linear behavior of Tc{gH) up to gH of the order 
of 1 GeV^, while at the moment there are not data for 
larger values of gH. Nevertheless the existing data are 
once again well interpolated by a linear fit. An extrapola- 
tion of the fit to larger values of the applied field strength 
implies that there exists a critical value of the applied 
field, gHc, such that the deconfinement temperature is 
zero, and for gH > gHc the system is deconfined at any 
temperature. Moreover, the Lattice results lead to the 
conclusion that the lowering of the deconfinement tem- 
perature is a peculiarity of the non abelian theories p^ . 
A possible explanation of these facts, as suggested by the 
authors in Refs. [H, is that the QCD vacuum be- 

haves as a relativistic color superconductor, whose color 
superconductivity is destroyed by a strong enough exter- 
nal magnetic field. 

In our calculations, we use the Nambu-Jona Lasinio 
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PNJL model in what fol- 

231, 123, M, M M M M M, 

as an effective model of QCD. 
The PNJL model has been widely used in the literature 
to depict several aspects of the QCD phase diagram. In 
the PNJL model, one has not dynamical gluons and the 
thermodynamics of the Polyakov loop L is driven by a 
temperature dependent effective potential U for L, the 
latter being coupled minimally to quarks via the QCD 
covariant derivative. Instead to use the usual hard cutoff 
regularization scheme, we consider here a non-local inter- 
action [H, m, im, S m, Si • we take an instanta- 
neous, rotationally and translationally invariant interac- 
tion, leaving the more general case of a non instantaneous 
interaction to a future study. The non local interaction 
is modelled in momentum space by a form factor which 
depends only on the magnitude of the 3— momentum and 
on a mass scale A. This approach has the advantage that 
it docs not artificially cutoff large momenta, thus making 
the model suitable for field strengths gH larger than the 
natural momentum scale A of the form factor. 

Since the PNJL model does not describe dynamical 
gluons it cannot be used to predict the response of the 
pure gauge system to an external field. Therefore in 
this paper we use the Lattice QCD result on the de- 
pendence of the dcconfincmcnt temperature of the pure 
gauge theory on oH as an input. This dependence has 
the form [iH, HI, |l3 



Tr 



= 1 — a- 



(1) 



where Tc is the deconfinement temperature at zero field 
and a is a dimcnsionless constant. However we arc able 
to compute the effect of dynamical quarks on deconfine- 
ment, and compare them with Lattice QCD results. Re- 
markably enough, our findings on deconfinement temper- 
ature with quarks are in agreement with Lattice QCD 
calculations. In particular, we will show that the linear 
dependence depicted in Eq. ([T]) is still valid in presence of 
quarks, as it is observed in Lattice measurements. Even 
if the larger part of our results are obtained in the physi- 
cal limit in which = 135 McV and for a specific choice 
of the form factor, we have checked the stability of our 
results by changing the analytical expression of the form 
factor in the chiral limit (finite quark masses are not ex- 
pected to play a relevant role in this context). 

The plan of the paper is as follows. In section II, we 



review the PNJL model with non local instantaneous in- 
teraction, and discuss the quark spectrum in the presence 
of an abelian external chromomagnetic field. In section 
III, we present our results on chiral and deconfinement 
transitions, as well as the phase diagram of the model in 
the gH — T plane and the equation of state of the PNJL 
matter in the external field. Finally, in section IV we 
draw our conclusions. 

II. THE PNJL MODEL IN EXTERNAL 
ABELIAN CHROMOMAGNETIC FIELD 

Our investigation concerns the response of the PNJL 
vacuum to an external abelian chromomagnetic field 
specified by the gauge potential 



An 



(2) 



corresponding to a field along the positive 2:-axes. The 
quark Lagrangian density is given by p^ . [2^ . [33| 



(3) 



In the above equation we have omitted the pure gauge 
contribution cx since it is a constant and thus not 
relevant for the dynamics; is the quark field with 
Dirac, color and flavor indices (implicitly summed), m 
corresponds to the bare quark mass matrix: we work 
with two flavors and we assume from the very begin- 
ning 771 u = rrid- The covariant derivative is defined as 



iAfj^ with 



gAlTa and r„ a = 1, 



being the SU{i) color generators with the normalization 
condition TY[Ta,Tp\ = Sab- 

In Eq. ^ L, L correspond to the normalized traced 
Polyakov loop and its hermitian conjugate respectively, 
L = TrW/Nc, L = TrV^V^c with 



\ j Aidr^ = exp (1/3^4) 



iAo 



(4) 

and 13 = l/T. In the above equation wc have implic- 
itly assumed that the thermodynamics of the Polyakov 
loop is described in terms of a constant and homogeneous 
background field Aq. Even if this choice leads to inter- 
esting agreement of the PNJL model with lattice calcu- 
lations, more sophisticated models with inhomogeneous 
background field have been investigated, see for exam- 
ple [H]. 



The term U[L, L, T] in Eq. ^ is the effective potential for the traced Polyakov loop; it is built by hand in order to 
reproduce the pure glue lattice data of QCD |19i.i2Q.i2ll.[26l? ]. In this paper, we work in the Polyakov gauge in which 

(5) 

with ^3, ^8 real parameters. Since we are interested to the case of zero baryon chemical potential we set (L) = (L) 
from now on. This choice implies = and thus we are left with only one parameter, £3 = £. Moreover, we adopt 
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the following logarithmic form [2l| 



U[L,L,T] = 



LL + b{T) log [1 - &LL + 4(L3 + L^) - iiLLf] 



with 



Numerical values of the coefficients are as follows [2l| : 

ao = 3.51 , ai = -2.47 , 03 = 15.2 , 63 = -1.75 



(6) 



(7) 



(8) 



r 



The parameter To in Eq. ([6]) sets the deconfinement scale 
in the pure gauge theory. Lattice results show that the 
deconfinement temperature of the pure glue system is a 
linear and decreasing function of \/gH [1^, , at least 
for small values of ^/gH/Tc, with Tc the deconfinement 
temperature at zero field. For ^fgH larger than a critical 
value of the order of 1 GeV, the critical temperature of 
the pure glue system as computed in Refs. [3, are 
consistent with zero. Inspired by these results, we make 
the following ansatz for Tq in Eq. ©: 



To 



270 MeV X 



1 - 



(9) 



with 6 denoting the unit step function, and we take 
\/gHc = 1.2 GeV according to the estimates of [Tsl. [l6j. 

Finally, £4 in Eq. ([3|) represents the lagrangian den- 
sity for the four fermion interaction. If we define Sa = 
f d^xLj as the interaction action then [2^. [38l. [39l. |40|. l4l]. 



Si^G j d^x [{q{x)q{x)f + {q{x)i^^Tq{x)f] , (10) 
with the dressed quark field defined as 

q{x) = f d^y F{x - y)^{y) , (11) 



where F{r) is a form factor whose Fourier transform /(p) 
satisfies the constraint /(p) ^ for p ^ 00, p being the 
3-momcntum. In this paper, we use the Lorcntzian form 
factor 



fip) 



1 



(12) 



In the above equation N ~ 10, p — \p\ and A ~ 
684.2 MeV. Moreover we use m = 4.46 MeV and G = 
2.33/A^ [2^. By these numerical values one has /^r ~ 
92.3 MeV and the pion mass = 135 MeV, as well as 
the zero momentum quark mass M„ = 335 MeV. The 
advantage of using a form factor in momentum space is 
that it does not lead to the introduction of an artificial 
momentum cutoff as in the usual NJL model calculations. 



In order to study chiral symmetry breaking we assume 
that in the ground state 



a = G{q{x)q{x)) ^ 



(13) 



where a summation over flavor and color is understood. 
In what follows we consider the system at finite tempera- 
ture T in the volume V. This implies that the space-time 
integral is / d'^x = Jq dr J d'^x with f3 = 1/T. In order 
to define the thermodynamical potential at 7^ 0, we 
observe that at iJ = the momentum space mean field 
PNJL action reads [13 



S = 



d^p 
d^p 



[4'{i^.p''-i^.Ani^] 

J{pf [2a ^{p)xP{p)\ 



(27r) 

2 

PV^ - PVU[L,L,T] , (14) 



with V denoting the quantization volume and A^^ = 
gA'^Ta- We introduce the mean field, momentum de- 
pendent, constituent quark mass M{p): 



M{p) EE m - 2af(p) = m + S(p) , 



(15) 



where S denotes the momentum dependent proper quark 
self energy. 

In order to compute the effective potential for a and L 
we need to know the quark spectrum. The computation 
of the latter in the case of non-local interaction is a non- 
trivial task, even in the mean field approximation: As a 
matter of fact, the Dirac equation for the quarks in the 
external field is an integro-diffcrcntial whose solution is 
beyond the scope of the present study. For this reason, 
in this paper we will make a simple ansatz for the quark 
spectrum of the non-local theory, arguing its specific form 
from that of the local NJL model. We recover the local 
NJL model by taking ^ 00 in Eq. (fT2|) . Hence, the 
larger the value A^, the better is our assumption. 
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The thermodynamical potential per unit volume at _ff = can be obtained by integration over the fermion fields 
in the partition function of the model: 

n = «|L.i.rl4-i-i:/^T,i„g5:^. (16, 

where the sum is over fermion Matsubara frequencies a;„ — T:T{2n + 1), and the trace is over Dirac, flavor and color 
indices. The inverse quark propagator is defined as 

S-\ioJn,p) = [{iujn+iA4)jo-l -P- M(p)](g)lf . (17) 

Performing the trace and the sum over Matsubara frequencies the effective potential for L, a at H ~ reads 

n = U[L,L,T] + ^~2N.NfJ^E 

-mfT J ^ log [1 + 3Le-f'^ + ZLe-^P^ + e'^/^^] , (18) 

where E = sjjp- + M(j)f- . 



The quark spectrum in absence of external field is given 

by 

< ^ ±i^±\/p2+M(p)2 , (19) 

= Tii^^v" + M{pY , (20) 

el = ^^TmW , (21) 

where the upper (lower) sign corresponds to particles (an- 
tiparticlcs). 

In presence of the external field ^ the red and green 
quark spectra is modified. If 54 = then it is well known 
that the effect of the applied chromomagnetic field on 
the quark spectra is both diamagnetic and paramagnetic. 
The former corresponds to the quantization of the quark 
motion in the plane orthogonal to the field direction. The 
latter differentiates the spin up and spin down motion 
by virtue of the usual term ^ a ■ H in the interaction 
hamiltonian. 

If 5*4 7^ but the interaction kernel is local [which cor- 
responds to the choice F{x — u)~ S'^{x — u) in Eq. (fTT|) ]. 
then the quark spectra in the case H ^ Q is 

E-^^^pI+pI^^ + A'P , E^,=^p^+AP , (22) 
where 

pl.s^gH (^n+^+-^^ , n = 0,l,2,... (23) 

and M denotes the constituent quark mass. In the case 
of the non-local interaction, we assume that the quark 
dispersion laws for red and green quarks are given by 

K = ^pI+p1.s + m{pI+pIJ^ , (24) 

while the spectrum of the blue quarks, which are not cou- 
pled to the external field, is still given by Eq. (PT|) . The 



assumption in Eq. ((M]) is partly justified by the observa- 
tion that 5*4 is translationally and rotationally invariant, 
and it does not contain derivative couplings, and at the 
mean field level, it simply amounts to the replacement 
of a constant constituent quark mass with a momentum 
dependent one, see Eq. The assumption ([M)) leads 

to the correct spectrum in the case ^4 can be treated as 
a perturbation, that is in the case of large gH with re- 
spect to a. This assumption has been used several times 
in the literature for studying chiral symmetry breaking 
in a strong magnetic field [4^, |4^ . 

In Equation ((24|) the subscripts r, b denote the quark 
color; s — ±1 distinguishes the quark spin projection 
along the applied field; finally n labels the Landau level. 
We have not explicitly written the energy of the green 
quarks since the former depends only on the absolute 
value of the charge, thus Eg = E^. Notice that we have 
introduced a spin index also for blue quarks: this is done 
for notational convenience. We have explicitly written 
the dependence on p± in Eq. ([M)) : in the following we 
will write M{p) for simplicity, leaving understood the de- 
pendence on the true quark momentum +p\ ■ Keeping 
into account that red and green quarks have respectively 
Ta color charges given by ±1/2, the phase space integral 
is 




for red and green quarks, and 

for blue quarks. Taking into account that in the Polyakov 
gauge the effect of the Polyakov loop is merely a shift of 
the quark chemical potentials to a complex value we can 
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write the quark poles as 



where the upper (lower) sign corresponds to particles (an- 
tiparticles). 



+ M(p)2 



(27) 
(28) 
(29) 



The thermodynamical potential at H ^ can be obtained from the expression a.t H = 0, see Eq. ([T8|), by replacing 
the quark energies and the phase space integral as follows: 



C +U[L, 



S = ±l 



(2^) 



3^b 



-2Nf J2 T 

s=±l 

-2Nf J2 T 



s=±l 



(27r)3 
(27r)3 



log 1 + iLe 



3Le 



log 



1 + e 



2NjJ2t 



s = ±l 



(27r)3 



log 



1 + e 



(30) 



In Eq. ([30|l C denotes an irrelevant additive constant 
that makes the thermodynamic potential equal to zero 
at T = 0, (7 = for each value of gH. This subtraction 
of the perturbative part of the thermodynamical poten- 
tial makes the latter finite at every temperature. Equa- 
tion (|5D|) is in agreement with that of Ref. [i^l in the 
limit of the pure NJL model, the latter being recovered 
in the limit L — > 1 after subtraction of U{L, L, T). 

The last line of Eq. (|30|1 corresponds to the thermal 
blue quark contribution subtracted of a fictitious contri- 
bution of the quarks of the same color in the external 
field. The same contribution is added to the second line. 
This procedure is applied for the phase space integral is 
different for red and green quarks on one hand and blue 
quarks on the other hand, thus their contribution to the 
free energy cannot be collected under an unique integral 
sign. This procedure has the advantage that it allows us 
to write the dependence of D, on L instead of £, the latter 
appearing explicitly in the quark dispersion laws. 

In order to avoid confusion with existing literature, we 
stress that in this paper we shall refer to a as the order 
parameter of chiral symmetry breaking and not to the 
chiral condensate. The identification of the two quanti- 
ties is correct in the case of a local interaction (apart an 
irrelevant multiplicative constant), but in the case of a 
non-local interaction they are no longer equal. This has 
been already discussed in a transparent way in Ref. [36j . 
The difference arises from the very definition of a in 
Eq. (|13p in terms of the dressed quark fields q. This is still 
more evident by inspection of the gap equation and the 
self consistent equation for the chiral condensate. The 
former is obtained by the requirement that 17 has to be 
minimized by the physical value of cr. At T = 0, this 



leads to the relation 
cr = -2GNcNf 



(27r)3-'^^^ y^TJTW ' 



(31) 



where we have taken the limit gH = (for gH ^ 
the same reasoning applies) and L = 1 (that is pure 
NJL model). On the other hand, the chiral condensate 
is defined in terms of the trace of the fermion propagator 
and of the bare quark fields, 



(i/jip) = 2{uu) = —i 



d^p 



TtS{p) (32) 



where the trace is over color, flavor and Dirac indices. 
The previous equation leads to 



{uu) — —6 



d^p M (p) 
(2^ VP^+Mp 



(33) 



In order to take into account only the non-perturbative 
contributions to the expectation value of the operator uu 
in the ground state, we have to average over the non per- 
turbative fluctuations only. To this end, the perturbative 
contribution to (uu) has to be subtracted from the pre- 
vious equation. Thus the non perturbative contribution 
to (uu), that has to be identified with the true chiral 
condensate, is given by 



{uu)np — —6 



d^p 
(27r)3 



M{p) 



(34) 

which shows that a ^ 2G{uu). We will make use of 
Eq. dSlD in section III.C. 
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Before closing this section we recall the definition of 
the dimensionless susceptibility matrix [l^, [23, [13 ■ The 
dimensionless curvature matrix of the free energy around 
its global minima C is given by 



C 



Cmm Cml 
Cml Cll 



with matrix elements defined as 



Cmm 
Cll 
Cml 



1 d^n 
TAdAP ' 
1 d^n 

1 d^n 
TA^d^dM 



(35) 

(36) 
(37) 
(38) 



and A'l = M{p — 0). The susceptibility matrix x is 
computed as the inverse of the curvature matrix C . We 
have 



X 



Xmm Xml 

XML XLL 



(39) 



Here xmm and xll denote respectively the dimension- 
less susceptibilities of the constituent quark mass and of 
the Polyakov loop. 



III. RESULTS 

We now discuss our results. In our model, we have in- 
troduced a form factor which mimics asymptotic freedom 
of QCD. Thus the hard momenta p ^ A are naturally 
cutoff in our theory. However from a numerical point of 
view we need to cut the integrations by hand at a certain 
value of momentum. In order to do this safely we ana- 
lyze the numerics of the form factor specified in Eq. (fT^ 
and notice that /(A) = 1/^2, /(2A) w 9.8 x lO"-* thus 
/(2A)//(A) w 1 X lO"'^. This means that our interaction 
is ineffective at momenta larger than 2A. Therefore, it is 
a good approximation to implement into the numerical 
integrals the condition p < 2A, which translates into 



< 4A^ 



(40) 



for red and green quarks. In order to satisfy this condi- 
tion, we require firstly 



p1.s < 



1 s 

- ± - < IntegcrPart 

2 2 - ^ 



4A- 



21 



MAX 



(41) 

where the upper (lower) sign corresponds to the case s = 
1 (s = —1)- The allowed values for n lie in the range 
{0, 1, . . . , Nmax - 1} for s = 1 and in {0, 1, ... , Nmax} 
for s = — 1. For gH > 4A^ we have Nmax ~ 0. For each 
allowed value of n we require secondly 



pi<9H 



,1s 

Nmax - I " + 2 ^ 2 



(42) 



In the following, we show results obtained with the choice 
given in Eq. ([1T|) . We have explicitly verified the stabil- 
ity of our numerical results by doubling the value of the 
numerical cutoff, thus using a maximum value of n which 
is four times larger than that in Eq. (|4T|) . 



Chromomagnetic catalysis of chiral 
condensation 



In Fig. m wc plot the zero momentum constituent 
quark mass Af„ = M{p = 0) as a function of the exter- 
nal field at zero temperature. The qualitative behavior 
of Mu is in agreement with that obtained previously for 
NJL in a rnagnetic field [i^ as well as in a chromomag- 
netic field |47l [49l . [sot with a proper time regularization. 
We notice that the net effect of gH is to increase the 
value of Af„ with respect to its value at gH = 0. This 
behavior is the so called (chromo)magnetic catalysis and 
it has been interpreted in terms of dimensional reduc- 
tion [H, [13, \M, m, MMEM- 

It is interesting to see how dimensional reduction man- 
ifests in the model at hand. For gH/ A^ ^ 1 the number 
of Landau levels that enter into the gap equation for cr, 
namely Nmax defined in Eq. (j41[) . is extremely large. As 
a consequence p^ can vary in a wide interval, see Eq. (|42p . 
In this case the motion is effectively three dimensional. 
As the magnitude of the external field is increased Nmax 
decreases. This is not an artifact of the model but a mere 
consequence of the existence of a form factor in the inter- 
action that naturally cuts off the hard momenta p ^ A. 
As Nmax decreases and eventually reaches the unity, 
the only possible value of p^. for spin up quarks is zero; 
regarding spin down quarks pi is zero for n = 1, while 
varies in the range (0, gH) for 71 = 0. Thus the motion in 
the lowest Landau level takes place in a two dimensional 
momentum space for spin up quarks, that is in the plane 
orthogonal to the applied field. The same is true for spin 
down quarks in the first excited Landau level. Finally 
the motion of spin down quarks in the lowest Landau 
level, with p\ = 0, takes place in an one dimensional 
momentum space. 



B. Deconfinement and chiral transitions 

In Fig.[2l we plot cr, normalized for each value of gH to 
its zero temperature value, the Polyakov loop and the two 
physically relevant susceptibilities, xmm and xll, as a 
function of temperature for three values oi gH. We iden- 
tify the deconfinement (chiral) crossover with the peak of 
the Polyakov loop (chiral) susceptibility xll (xmm)- An 
interesting feature of the model at hand is that the decon- 
finememt crossover becomes a first order transition when 
the strength of the applied field is larger than a critical 
value. This is clear both from the behavior of (L) against 
T and from that of xll- As a matter of fact, at gH = 
the Polyakov loop increases smoothly as the temperature 
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is increased. At large values of gH the Polyakov loop is 
consistent with zero at small temperature, and becomes 
suddenly non zero at large temperature. The value of the 
critical temperature depends on gH . Moreover, xll has 
a broad and not pronounced peak at gH = 0, signaling 
the transition is actually a smooth crossover (in Fig. the 
value of Xll at gH = is multiplied by a factor of 10 
for a better comparison with the other data). As gH is 
increased, xll develops a pronounced peak, the larger 
the value of gH the larger the height of the peak. Thus 
the crossover becomes a true first order transition. 

On the other hand, from the data we obtain for xmm, 
we infer that the chiral transition remains a crossover 
at each value of gH we have considered. The effect of 
gH on chiral symmetry breaking is twofold. At zero 
temperature we find that the larger the magnitude of 
gH the larger the magnitude of the chiral condensate. 
This is in agreement with the aforementioned scenarios 
of chromomagnetic catalysis. At finite temperature the 
scenario changes depending on the value of gH. The chi- 
ral crossover temperature slightly decreases at small gH, 
but increases at large gH . 

In Fig. [31 we plot the critical temperature for de- 
confinement as a function of the applied field strength 
\fgH . We identify the deconfinement temperature with 
the peak of the Polyakov loop susceptibility. In the fig- 
ure, the black rectangles are the results of our calcula- 
tions. Interestingly enough, we can fit the data with a 
linear function of \J gH . We define Tc as the deconfine- 
ment temperature at gH = with dynamical quarks, 
namely Tc ~ 217.6 McV. Then we find the best fit to our 
data as 

^ = 1-0.176^, (43) 

with a linear regression coefficient i?^ = 0.996. Eq. ^ 
can be considered the main result of our investigation. 
The linear dependence of the deconfinement tempera- 
ture on the square root of the external field strength has 
been noticed for the first time in [l^, [l^ within lattice 
calculations in the pure gauge theory. Existing Lattice 
data sustain this picture even in presence of dynamical 
quarks [l3| , at least for \/ qH of the order of 1 GeV (see 
for example Fig. 7 of Ref. [13] )■ In our model we use the 
linear fit of [TsTIl^ on the deconfinement temperature of 
the pure gauge as an input. On the other hand, the linear 
dependence of Tc{gH) on ^/ gH with dynamical quarks is 
an output of our calculations. Hence our result strength- 
ens the idea that the PNJL model captures the essential 
characters of the deconfinement transition. 

It is also interesting to notice that there exists a critical 
field, that we denote by \/gHc, at which the critical tem- 
perature for deconfinement vanishes. This is a mere con- 
sequence of the linear dependence depicted in Eq. (j43p . 
By requiring TdgHc) = 0, wc find from Eq. (|43[) 



\/gHc — 1.30 GeV , with dynamical quarks , (44) 



in good agreement with the value \J gHc = 1.6 GeV ob- 
tained by an extrapolation of the Lattice data [l3| . The 
existence of a critical chromomagnetic field for deconfine- 
ment is supported by the results of Ref. [i^ : As a matter 
of fact, it has been shown there that for strong enough 
chromomagnetic fields the color superconductivity, thus 
a deconfined phase, occurs even at zero chemical poten- 
tial. 



C. Changing the form factor in the chiral limit 

In this section, we study the effect of choosing dif- 
ferent form factors on the results discussed above. To 
begin with, we show that the oscillations of the function 
Mu (gH) in Fig. [1] is definitely due to the sharpness of 
the form factor: choosing a smoother form factor the os- 
cillation disappears. For our purposes, it is enough to 
take the chiral limit and to focus on the lorentzian form 
factors /jv defined by the equation 

fN{p) = ^==1=== , P^\p\ ■ (45) 

^1 + (p2/A2)JV 

In the previous sections, we have shown the results cor- 
responding to = 10. We now consider other cases and 
compare them to the case A^ = 10. As A^ is decreased, the 
form factor becomes a smoother function of p/A. This 
is shown in Fig. |4] where the form factors corresponding 
to A^ = 10,3 and 1.5 are shown for a common value of 
A = 700 MeV. 

In order to properly fix the parameters G and A for 
any given value of A^ in Eq. (j45p . we require that the 
model reproduces the experimental value of /,r = 92.3 
MeV and a phenomcnological value of the constituent 
quark mass at zero momentum and at zero field. We 
take the latter in the range 300 — 400 MeV, correspond- 
ing to an expectation value of |(t| in the range 150 — 200 
MeV (for recent estimates of this quantity see for exam- 
ple [Sy, [131). Once the parameters are fixed at a given 
value of A", we measure the value of the subtracted chiral 
condensate defined by Eq. ((M]) . taking care it is in the 

1 /3 

allowed phenomcnological range 190 MeV < — (uu)^p 

< 260 MeV [53. This constraint on {uu)](p can be ful- 
filled only for relatively large values of A^, namely N > 5, 
within the instantaneous form factor approach [i^, [3] ; 
on the other hand, more refined calculations based on 
non-instantaneous form factors can fulfill the above con- 
straint [H, [3| also for smoother form factors. Therefore 
the picture we draw in this section has to be considered 
purely qualitative, at least for relatively small values of 
the parameter A^. The numerical values of the parame- 
ters are collected in Table [Ij 

In the upper right panel of Fig. |4l we plot the quark 
self energy at zero momentum as a function of the applied 
field strength. Wc have shown results for A^ = 10 (solid 
line), A^ = 3 (dashed line) and N = 3/2 (dot-dashed 
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(J (MeV) 


A (MeV) 


GA^ 


1/3 

-{uu)j!jp 


fio 


-160.2 


660 


2.36 


256 


h 


-160.2 


679 


2.49 


269 


h 


-170.2 


674 


2.62 


282 


fz/2 


-200 


614 


2.51 


312 



TABLE I: Parameters for the various form factors examined. 
For all the values of A'^ the inputs are /,r = 92.3 MeV and the 
value of o; A, GA'^ and —{uu)\li^ are outputs. 

line). As anticipated the effect of a smoother form fac- 
tor is to reduce (and eventually definitively damp) the 
oscillations observed for the case N = 10. 

We have checked the stability of our main results on the 
deconfinement temperature in external field by changing 
the form factor to a smoother one in the chiral limit. 
The results are summarized in the lower panels of Fig. HI 
where we plot the expectation value of the Polyakov loop 
computed for the form factor fsip), and the deconfine- 
ment line. In this case the deconfinement temperature at 
zero field is = 213 MeV. The dashed line corresponds 
to the linear fit to the data, specified by the equation 

^-1-0.161^, (46) 

with a linear regression coefficient i?^ = 0.99. In this 
case, the critical field corresponding to a vanishing de- 
confinement temperature is 

V5^^= 1.31 GeV , (47) 

in agreement with our previous estimate obtained us- 
ing TV = 10 and with quarks in the physical limit, see 
Eq. (glD. 

D. The phase diagram in the gH — T plane 

We summarize our results in a phase diagram in the 
gH — T plane, see Fig [5j In the figure, the dashed line 
denotes the deconfinement crossover, which becomes a 
first order transition as gH is larger than a critical value, 
see the bold line in Fig. [S] Both the dashed and the solid 
lines are the results of the fit given in Eq. (|43)) and are 
identified by the position of the peaks of the Polyakov 
loop susceptibility. They are an output of our calcula- 
tions once the input Eq. ^ is given. Our results show 
that there exists a critical endpoint in the phase diagram, 
with coordinates given by 



(^,f) = (3.1,0.5), (48) 

that is, 

(x/gHE, Te) = (0.653, 0.108) GeV . (49) 

The critical end point coordinates have been computed 
by the investigation of the Polyakov loop temperature 



dependence. At a given value of gH, we compute L as a 
function of temperature. If the function L{T) is contin- 
uous then we identify a crossover. On the other hand, if 
L{T) shows a discontinuity, then we identify the critical 
temperature by the temperature at which the discontinu- 
ity occurs, and identify the transition with a first order 
transition. By iteration of this process for several values 
of gH, we are able to determine the critical end point 
coordinates. The uncertainty of our final result depends 
only on the size of the interval of \fgH , around gHE , in 
which we look for L{T), and we can estimate it to be of 
the order of 10 MeV, both on T and on ^/gH 

In Fig. [5] we draw also the chiral crossover line. It is 
obtained by the position of the peaks in the chiral sus- 
ceptibility. For gH < gHE, the effect of the applied 
field is to slightly lower the transition temperature. For 
gH > gHE, the slope of the chiral crossover line changes 
sign thus opening a large window with dcconfincd mat- 
ter but with chiral symmetry spontaneously broken. The 
existence of the valley in the chiral critical line around 
gHE could be an artifact of the specific functional choice 
of the form factor. On the other hand, the qualitative 
behavior we have depicted should be quite robust since 
the dimensional reduction and the consequent chromo- 
magnctic catalysis exist in the model at hand (see the 
discussion at the beginning of this section) independently 
on which form factor is chosen. 



E. Equation of state in the external field 

For completeness, we compute in this section the equa- 
tion of state p = p{£) of the PNJL model for different 
values of the strength of the external field. This kind of 
computation is interesting because the equation of state 
is accessible to Lattice measurements, hence our results 
can be directly checked. The various thermodynamical 
quantities are defined as follow: 



p = 


-n , 


(50) 


s = 




(51) 


e = 


-p + Ts , 


(52) 


Cv = 




(53) 


cl - 


dp s 
de Cv 


(54) 



In the above equations p, s, e, Cv and denote the 
pressure, the entropy, the internal energy density, the 
specific heat and the squared sound velocity, respectively. 
We normalize the pressure with the convention p{T = 
0) = for each value of gH. 

The results of our numerical calculations are shown in 
Figs. [6] and [7] where we plot the pressure and the energy 
density, the specific heat, the squared sound velocity and 
the equation of state for gH = and for gH = 0.8 GeV^. 
From the qualitative point of view the thermodynamical 
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quantities behave like the gH = case for gH < gHE, 
and Uke in the case gH = 0.8 GeV^ for gH > gH^- For 
this reason, it suffices to plot the various quantities only 
for these two values of the applied field strength. 

IV. CONCLUSIONS 

In this paper, we have investigated the response of the 
PNJL vacuum to an external abelian chromomagnctic 
field. Our interaction is non-local but instantaneous, 
with a momentum-dependent Lorentzian form factor. We 
have used as an input the Lattice-deduced dependence of 
the deconfinement temperature of pure gauge theory on 
the strength of the applied field gH, and we have com- 
puted the effect of dynamical quarks on deconfinement. 

We can summarize our results as follows: An external 
abelian chromomagnctic field inhibits color confinement 
even the in presence of dynamical quarks. Moreover the 
dependence of the deconfinement temperature on \/gH is 
linear. We have compared our results with existing Lat- 
tice data, the latter exploring values gH up to the order 
1 GeV^ in the case of full QCD, finding qualitative agree- 
ment. An extrapolation of our data at larger values of gH 
leads to the existence of a critical field gHc above which 
quark matter is deconfined even at small temperatures. 
We have checked the stability of this result by changing 
the analytical form of the momentum-space form factor. 
We have drawn the phase diagram of the model in the 
gH — T plane in the physical limit. Finally we have com- 
puted the equation of state of the PNJL matter in the 
external field. 

The results on the deconfinement temperature are en- 
couraging: They suggest that the PNJL model captures 
the essential physics of deconfinement. On the other 
hand we find that the external field favors spontaneous 
chiral symmetry breaking, as in the NJL model (chromo- 
magnctic catalysis). The latter result is in disagreement 
with Lattice, where a single transition temperature is 



observed at which both the Polyakov loop and the chiral 
susceptibilities show pronounced peaks. Thus, Lattice 
supports a scenario in which deconfinement and chiral 
restoration occur at the same temperature for any value 
of gH] as a consequence spontaneous chiral symmetry 
breaking should be inhibited by the external field and 
not catalyzed. Probably, the discrepancy arises mainly 
from the fact that our interaction is instantaneous, that 
is it depends only on the three-momentum. 

A comparison with QCD in st rong magnetic fields can 
be helpful. The authors of Ref. [4^ have shown, within 
an improved ladder approximation, that in QCD a strong 
magnetic field, B, leads to a dynamical mass mq{B) 
which is lower than mq{B = 0), if eB is in the range 
^QCD ^ < (10 TeV)2. This is mainly due to the 
running of Ug at the scale eB, 



Moreover, in the momentum region relevant for dynami- 
cal symmetry breaking, that is — < eB, the 
magnetic field acts as a cutoff in the whole 4-momcntum 
space; on the other hand we have assumed that all the 
energy region is relevant for dynamics (we have no cutoff 
on energy), and our coupling constant G is fixed, once 
for all, at the value corresponding to zero field. These 
two effects combined together probably lead to an over- 
estimate of the chiral condensate, as argued in Ref. [4^. 
Investigations based on non-instantaneous interaction as 
well as on running coupling are now under investigation 
and will be discussed in a forthcoming paper. 
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FIG. 1: Zero momentum quark mass as a function of the applied chromomagnetic field at zero temperature. 
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FIG. 2: Expectation values of a, normalized to its zero temperature value, and of the Polyakov loop as a function of temperature 
(in units of the zero field deconfinement temperature Tc = 217.6 MeV), for three different values of gH. In the lower panel the 
chiral and the Polyakov loop susceptibilities are shown. 
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FIG. 3: Deconfinement temperature against \/gH for the model with dynamical quarks. Tc is defined as the deconfinement 
temperature at zero field, Tc = 217.6 MeV. The black rectangles correspond to the results of our calculation; the dashed line 
is a linear fit to the data. 




T/Tc T/Tc 

FIG. 4: Upper left panel: form factors for three different values of A'^ for a common value of A = 700 MeV. Upper right panel; 
quark self energy at zero 3-momentum as a function of the applied field strength. Lower left panel: expectation value of the 
Polyakov loop computer for three different values of the applied field strength, in the chiral limit and for the f^ip) form factor. 
Lower right panel: deconfinement critical line computed in the chiral limit and for the faip) form factor. In the lower panels 
the deconfinement temperature at zero field is Tc = 213 MeV. 
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FIG. 5: Phase diagram of the PNJL model in the gH — T plane. Mass scale is given by Tc = 217.6 MeV corresponding to the 
deconfinement temperature of the model with dynamical quarks and at gH = 0. Dashed line corresponds to the deconfinement 
crossover; solid line denotes the first order deconfinement transition. Dot-dashed line corresponds to the chiral crosser. 
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FIG. 6; Upper left panel: pressure and energy density, normalized to T*, against temperature in the case gH — 0. Upper 
right panel: pressure and energy density, normalized to T*, against temperature in the case gH = 0.8 GeV^. Lower left panel: 
dimensionless specific heat Cv /T^ as a function of temperature for gH = and gH = 0.8 GeV'^. Lower right panel: squared 
sound velocity against temperature. In the figure the deconfinement temperature at zero field is Tc = 217.6 MeV. 



